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Abstract. - Experimental results on amorphous ZrAlCu thin film growth and the dynamics 
of the surface morphology as predicted from a minimal nonlinear stochastic deposition equation 

S are analysed and compared. Key points of this study are (i) an estimation procedure for 
coefficients entering into the growth equation and (ii) a detailed analysis and interpretation of 
the time evolution of the correlation length and the surface roughness. The results corroborate 
the usefulness of the deposition equation as a tool for studying amorphous growth processes. 

P 
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' Introduction. - During the last decade the study of the kinetics of surface growth pro- 



cesses has attracted considerable interest (cf. the reviews ) . The dynamics of the surface 
morphology e.g. in amorphous thin film growth is dominated by the interplay of roughening, 
smoothening, and pattern forming processes. On the microscopic level, these processes are 
governed by the highly complex and only partly understood interaction of the depositing 
particles with the already condensed surface atoms. Despite the complexity of the growth 
processes on the atomic scale, experiments on the slightly coarser mesoscopic scale typically 
reveal some sort of regularity of the surface morphology with some superimposed small-scale 
stochastics [|[ [|. This, in turn, indicates that the machinery of coarse-grained continuum 
models based on phenomenologically motivated stochastic growth equations Q is a useful tool 
for the understanding and interpretation of the growth dynamics. In particular, amorphous 
thin film growth represents an attractive testing ground for the validation of such phenomeno- 
logical models; this is mainly due to the spatially isotropic nature of the amorphous structure 
at this scale and the lack of long range ordering. 

Our objective is a detailed comparison between a stochastic nonlinear evolution equation 
for amorphous thin film growth p| and experimental results on the surface morphology of 
ZrAlCu films prepared by physical vapor deposition that are analysed using scanning tunneling 
microscopy (STM). Using the aforementioned theoretical approach, we develop a method to 
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estimate the phenomenological parameters that is based on the evolution for short times. We 
also investigate whether the finite apex angle of the cone-like STM tip used in the experiment 
does affect the comparison between theory and experimental findings. In particular, we 
demonstrate that the results for both the surface roughness w (i.e. the root mean square 
deviation of the relative height fluctuations) and correlation length R c (i.e. the typical length 
scale over which height fluctuations are correlated) of our full nonlinear growth equation yield 
very good agreement with the presented experimental data, with only minor modifications due 
to the finite apex angle. This clearly corroborates the usefulness of our modeling approach. 

Experiments. - The glassy ZrAlCu films (composition ZrgsAly.sC^r.s) are prepared in ultra 
high vacuum by physical vapor deposition on oxidized Si wafers using a total deposition rate 
of F = 0.79nm/s (electron beam evaporation of the pure elements, each source independently 
rate controlled). Due to the geometrical arrangements, the particle flux is almost perpendicular 
to the substrate which is rotated during deposition (for a sketch of the experiment setup 
cf. fig. Q). The surface profiles of the films (sample area (200 x 200)nm 2 ) are analysed 
in situ by scanning tunneling microscopy using a tungsten tip with an apex angle of ca. 
40° (for further details cf. Rcf. Q). The film composition is verified by Auger electron 
spectroscopy. Additionally performed X-ray diffraction and differential scanning calorimetry 
(DSC) measurements ensure the amorphicity of the investigated films. From the STM data, 
the correlation length R c and the surface roughness w (for their definition, see below) are 
determined for various layer thicknesses up to 480nm. The experimental results correspond to 
the diamond symbols depicted in fig. || and || below. At large layer thicknesses (> lOOnm) 
each such symbol represents a measurement with a different sample. 

Theoretical modeling. - To accomplish a theoretical description of the time-evolution of 
an initially flat surface H(x,t), where H denotes the z-coordinate of the growing surface at 
the position x = (x,y) and time t (cf. also fig. [l]), we take advantage of the well-established 
phenomenological approach that is based on stochastic nonlinear partial differential equations 
|, i-e. 

d t H = G{H,VH) + F + v . (I) 

In eq. (|fj), G denotes a functional of the surface height and its local derivatives. The detailed 
functional form of G depends crucially on the considered experimental setup and the kinetics of 
the deposition process. F denotes the mean deposition rate and rj(x, t) is the related deposition 
noise that determines the fluctuations of the deposition process around its mean F . These 
fluctuations are assumed to be Gaussian white, 

(r 1 (x,t)} = 0; (r 1 (x,t)r ] (y,t')) = 2D6 2 (x-y)5(t-t') (2) 

where the brackets denote ensemble averaging. It proofs useful to introduce the height profile 
h(x, t) — H(x, t) — Ft. If the deposition process has any excess velocity, there is a nonlinear 
relation between the mean growth or layer thickness (H) and time, i.e. (H) (t) — Ft + (h) (t) 
with (h) (i) ^ 0. 

The simplest nonlinear amorphous thin film growth equation of the functional form ([!]) that 
incorporates (i) the physical symmetries such as rotation and reflection invariancc in the plane 
perpendicular to the growth direction, cf. fig. |l|, (ii) no particle desorption, and (iii) the 
potentiality of local density variations in the amorphously grown material yields in terms of a 
low-order expansion of G in the gradients of the surface profile h(x, t) the result 

d t h = ai \7 2 h + a 2 V 4 h + a 3 V 2 (Vfr) 2 + a 4 (Vfr) 2 + a 5 det ( ® x £ ^Ja? ) + V- ( 3 ) 
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with di(i = 1,...,5) being scalar material-dependent coefficients. In eq. (||), the first and 
the fifth term on the r.h.s. are directly related to the deflection of the initially perpendicular 
incident particles due to the interatomic forces between the surface atoms and the incident 
particles. For an indication of the relevance of this effect we refer to the recent experimental 
study in The coefficients ai and a 5 are determined by the relations a\ = —Fb and 
a 5 = Fb 2 , where b is the difference between the typical range of the interatomic forces and the 
equilibrium distance of the adatoms to the surface. Because b is small for the aforementioned 
experimental setup (typically of the order 10^ 1 nm) the term proportional to 05 in eq. (|^) 
can safely be neglected. The second and the third term on the r.h.s. of eq. (||) are related 
to the known microscopic mechanisms of (i) the surface diffusion suggested by Mullins [^| 
and (ii) the equilibration of the inhomogeneous concentration of the diffusing particles on 
the surface, as suggested in R Moreover, the coefficients and 03 are negative. The 
coefficient 03 reads 03 = — FZ /8, where I 2 is the mean square of the diffusion length of the 
particles. This characteristic dependence for amorphous growth is similar in nature to a term 
appearing for crystalline growth ||; it differs, however, due to the absence of an additional 
length scale, the typical height of crystalline terraces. The fourth term on the r.h.s. of eq. (|^) 
is of Kardar-Parisi-Zhang (KPZ) form (To). It is due to the potential dependence of the local 
density on the surface slope V/i, i.e. [p(V/i)] _1 = p Q 1 [l + (o 4 /F)(V/i) 2 + C((V/i) 4 )], with o 4 
being necessarily positive because of the additional volume increase caused by oblique particle 
incidence. Note also that a finite 04 results in a finite excess velocity. 

The experimentally detected correlation length R c (t) and surface roughness w(t) are deter- 
mined by the height-height correlation function 



C(r, t)= / i j d 2 x (h{x, t) - h) (h(x + f,t)-h)\ 



(4) 

\r\—r 



where h is the spatial average of the surface profile h = J d 2 yh{y 1 1) and L 2 the sample area. 
Specifically, R c (t) is given by the first maximum of C{r, t) occuring at nonzero r and the square 
of the surface roughness results by taking the limit r = in C(r, t), i.e. w 2 {t) — C(0,t). To 
integrate numerically the growth equation, a forward-backward finite difference method on a 
quadratic lattice combined with an Euler algorithm in time and periodic boundary conditions 
on a quadratic area [0, L} 2 have been invoked 

Linear model equation. - In order to make contact with the aforementioned experiment, we 
first extrapolate the values of the parameters a±, 02, and D. The initial stages of the deposition 
process (corresponding to layer thicknesses (H) < 240nm in the experiment), which started 
from a flat substrate, are dominated by the temporal evolution of the linear limit of eq. (^) , 

d t h = a 1 V 2 h + a 2 V 4 'h + r). (5) 

The growth of the Fourier modes of the surface profile is determined by dth(k, t) — a(k)h(k, t) + 
f]{k, t) where cr(fc) = — ai/c 2 +02fc 4 denotes the growth rate of the Fourier modes. Therefore, the 
wave number k c belonging to the maximum of the growth coefficient cr(fc) reads k c = \J a\/2a2- 
The correlation length R c (t) of the surface profile h(x,t) arising from the linearized equation 
(||) follows initially a i 1 / 4 -law, whereas it saturates at later stages into R c (t) — 7.0156/fc c = 
7.0156^/202/01 when the critical mode with the wave number k c dominates. Since a saturation 
of R c {t) with increasing layer thickness is also observed in the experiment for layer thicknesses 
larger than 200nm (cf. the diamond symbols in fig. |2|a), the ratio 02/01 can be roughly 
estimated. The surface roughness w(t) at later stages follows an exp[a(k c )t] = exp(—a 2 t/4a2)- 
behaviour when the critical mode dominates. By a comparison with the experimentally 
observed increase of w(t) (cf. the diamond symbols in fig. |^b) during the time interval when 
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the layer thickness (H) is between 30nm and 240nm the ratio a\/a 2 can be estimated. From 
the ratios 0,2/0,1 and a\/a2 the coefficients a\ and 02 can be deduced. The height-height 
correlation function C (r, t) that arises from the linearized equation ([5]) is determined by 



(2^) 2 J l K ' - ai k 2 +a 2 k 4 



(6) 

\r\—r 



From C(r, t) the surface roughness w(t) and the correlation length R c (t) can be determined. 
This makes it possible to determine more precisely the coefficients of eq. (0). The physical 

_2 _ -4 _ 2-2 — - 

dimensions of these coefficients are [ai] = I /t, [02] = I ft, and [D\ = h I ft where I denotes 
a length unit, t a time unit, and h a height unit, respectively. Therefore, ya^/ai is a length 
constant, \aifd\ | is a time constant, and yj Dj\a\\ is a height constant. Changing 02/0,1 by 
an arbitrary factor would change all lengths by the same factor. So we calculated w(t) and 
Rc(t) by means of eq. (^) with the approximately determined coefficients a\ and a 2 and an 
arbitrary D. Then, we changed R c , time t, and w by independent factors and \fa^fa\, |a2/a 2 |, 
and y/D/\a-\\ by the same factors, respectively, until R c (t) and w(t) were in accordance with 
the experimental result. By means of y/ 0,2/0,1, {0,2/ a1\, and y/D/\ai we evaluated ax, 02, and 
D. The parameters that fit the experimental results best, are found to read 

ai = -0.0826nm 2 /s, a 2 = -0.319nm 4 /s, D = 0.0l74nm i /s. (7) 

The dotted lines in fig. ||a and ||b depict the corresponding numerical simulations of eq. 

There is obviously good agreement with the experimental results for layer thicknesses 
(H) < 240nm. 

Next, we investigate possible effects caused by the STM tip. Although the detailed shape 
of the STM tip is known only roughly a reasonable model, which refers to the worst case 
situation, consists of a combination of a taper shank with an apex angle of ca. 40° and a top 
cone having a considerably wider apex angle of approximately a = 1 00°. Only the latter part 
is of relevance for the detection of the surface profile. The effect of the STM tip is that it maps 
the height profile h(x,t) from eq. (|^) on the "scanned" height profile h(x,t). The solid line in 
fig. |^b depicts the corresponding finding for this worst case of a STM tip with an apex angle 
of 100°. The major effect of the finite tip angle results in a slow down of the growth of the 
surface roughness w(t). This is caused by the scanning procedure that can not fully resolve 
the grooves of the height profile h(x, t). As can be deduced from fig. ^b, however, the linear 
growth eq. (|5|) (even including the tip angle) is not sufficient to describe the later stages of 
amorphous thin film growth. 

Nonlinear model equation. - We next investigate the physical relevance of the nonlinear 
terms a3V 2 (V/i) 2 and a^(VK) 2 on the growth dynamics, i.e. we analyse the nonlinear growth 
model 

d t h = ai V 2 h + a 2 V A h + a 3 V 2 (V/i) 2 + a 4 (V/i) 2 + 77. (8) 

The parameter choice that fits the experimental results best, is found to read 

a 3 = — O.lOnm 3 / s and 04 = 0.055nm/s, (9) 

where the parameters given in eq. (^) were used. The corresponding results for the correlation 
length R c (t) and the surface roughness w(t) of eq. are depicted in fig. || with (solid lines) 
and without (dotted lines) the inclusion of the effect of the STM tip. We now obtain good 
agreement with the experimental data also at larger layer thicknesses. Moreover, the finite 
apex angle of the tip yields only a minor effect (up to the largest experimental layer thickness 
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« 480nm), see the barely visible dotted line in fig. gb. As a general consequence, the two 
nonlinear terms result in a drastic slow down of the increase of the surface roughness w(t) 
with time or layer thickness. The larger the absolute values of 03 and 04 are, the stronger 
is this slow down. Nevertheless, the inclusion of both nonlincarities is necessary to obtain a 
consistent agreement with the experimental data. 

In particular, by setting a 4 = we need to decrease further a 3 in order to fit reasonably 
well the data. Nevertheless, the surface roughness w(t) is now either too large at large 
thickness (H) = 480nm, or too small at (H) = 360nm (not shown). The inclusion of 
the finite tip geometry does not cure this finding. Omitting 03 we find an overshoot for 
the roughness w(t) above (H) w 240nm. In addition, we note that the correlation length 
R c (t) ceases to exist above (H) w 300nm, because the first maximum of the hcight-hcight- 
correlation C(r, t) vanishes. These features are depicted by the dashed lines in fig. ||a and ||b. 
Fig. ||b also depicts that the experimental data for w(t) are strongly scattered at largest layer 
thickness. Thus, a precise estimation of the nonlinear coefficients is difficult. We estimate 
that the coefficient 03 ranges at most between — 0.13nm 3 / s < 0,3 < — 0.08nm 3 /s and likewise 

0. 04rtm/s < 04 < 0.07nm/ s. Consistent with the fitting procedure in fig. ^a and ||b we find 
that the crossover from linear to nonlinear growth behaviour sets in at (H) ps 240rtm. An 
inclusion of the small correction proportional to 0,5, see eq. (|^), does not quantitatively impact 
our results: For the largest layer thickness we find an improvement of maximal 1.0% for the 
roughness w(t) and 2.5% for the correlation length R c (t). 

Discussion. - The aforementioned extrapolation of the parameters a%, 02, 03, 04, and D 
also allows for additional microscopic estimates, (i) Since a\ = —Fb, the typical range b of 
the interaction between the surface atoms and the particles to be deposited, is about O.lnm, 

1. e. it approximately equals the size of the radii (0.2nm) of the surface atoms, (ii) Since 
03 = —Fl 2 /8, the diffusion length / must be in the range of l.Onm. This substantiates that 
the deposited particles experience a surface diffusion on a nanometer scale and do not just 
stick at the places where they hit the surface, (iii) If the particles arrive independently on the 
surface, the deposition noise is related to the particle volume fl and the mean deposition rate 
F by 2D = Fil, yielding (by use of eq. (0)) fl = 0.04nm 3 . This agrees within a factor of 
two with the averaged particle volume of ZrAlCu. (iv) The necessity of the inclusion of the 
term proportional to indicates that the local density of the growing film varies with the 
surface slope: On an inclined surface area the local density is decreased by p{Vh) = po/~f with 
7 = 1 + (a4/F)(Vh) 2 , where a\fF is in the range of 0.07. Note that these predicted finite 
density variations are physically compatible with the small diffusion length / of two to three 
atom diameters. At largest layer thickness (~ 480nm) this local density reduction (averaged 

/ — 1/2 

over the surface) possesses a mean 7 = 1.021 and a standard deviation I (7 — 7) J = 0.017. 

The maximum of 7 on the same surface (at L = 200nm) is typically of the order of 1.12. At 
the same time, 7 averaged over the whole film equals (H) /(Ft) = 1.01. 

Conclusions.- Starting from experimental data for amorphous ZrAlCu thin film growth, the 
phenomenological growth equation in ([}]), and using a step- by-step parameter identification 
procedure, we have shown that a quantitative agreement between theory and experiment is 
achieved. Based on our comparison between theory and experiment, we conclude that (i) eq. 

constitutes a valid theoretical model for amorphous thin film growth (at least up to the 
considered layer thicknesses) and (ii) that an interpretation of the data necessarily requires the 
inclusion of the two nonlinear terms V 2 (V/i) 2 and (Vft.) 2 in the stochastic growth equation (||). 
Several questions, however, remain open for future studies: (a) From the available experimental 
data it is not clear whether the correlation length R c (t) and the surface roughness w(t) saturate 
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for larger layer thicknesses, (b) The validity of the proposed approach needs to be tested more 
thoroughly by direct comparison with experimental data for spatio-temporal quantities such 
as the height-height correlation function, (c) A further experimental challenge presents the 
validation of the theoretically predicted local density variations. 
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Fig. 1. - Sketch of a physical vapor deposition experiment at normal incidence for amorphous film 
growth on a substrate. 
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Fig. 2. - The correlation length R c 
layer thickness (H), see below eq. 
parameters given in eq. (W). In fig. 



md surface roughness w (dotted lines), defined belo-vv eq. (jfy, vs. 
2|), are calculated from the linear growth equation (|5|) using the 
Jb we depict the influence of the STM tip with an apex angle of 



100° (solid line). The diamond symbols represent the corresponding experimental results. 




Fig. 3. - Correlation length R c and surface roughness w calculated from the nonlinear growth equation 
(@) using the parameters 03 = — O.Wnm 3 / s and 0,4 = 0.055nm/s with (solid lines) and without (dotted 
lines) consideration of the STM tip (with an apex angle of 100°). The parameters of the linear parts 
of the equation, 01, 02, and D, are given in eq. (Q). To demonstrate the physical relevance of the 
nonlinear growth law cc V 2 (V/t) 2 we depict for comparison the prediction that results by setting 
03 = (dashed lines without consideration of the STM tip). The diamond symbols represent the 
corresponding experimental results. 



